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Transverse Beam Dynamics including Aberration 
Effects in the Thermal Wave Model using a 
Functional Method 
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c/2 • We studied the transverse beam dynamics including aberration effects of sextupole and 
. ■ octupole perturbations in a thermal wave model. A functional integration method was 
used to calculate the first-order perturbation effects. We found that the model successfully 
explains a PARMILA simulation results for proton beams without space-charge effects in 
, CU a FODO lattice. 
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The thermal wave model (TWM) can be used to study the beam dynamics for rela- 
tivists charged particles pQ. This model was used to estimate the luminosity for a beam 
perturbed by the spherical aberration in a linear collider [2] and explain the halo forma- 
tion by using a Gaussian slit [3]. The transverse beam dynamics of charged particles is a 
promising branch of the TWM. In ref. [I], authors showed that the particle distribution 
was successfully described by the TWM in a thin lens approximation of a quadrupole mag- 
net with aberration. It can be extended to a general quadrupole lattice for a Gaussian 
beam by using a functional integration method if the space charge effects can be neglected 
[5]. However the work was restricted to a perfect quadrupole magnet. In this paper, we 
extended the TWM formulation for a transverse beam dynamics to include a sextupole 
and an octupole perturbations by using the first-order perturbation theory of a functional 
method. The model calculation was compared to the PARMILA [6] simulation result in a 
FODO lattice. We found that the model can explain the PARMILA results successfully if 
the perturbation effect is not too large. 

In TWM, the beam wave function of the relativistic charged particle can be described 
by the Schrodinger-type equation as following, 

d ip(x, z\ € 2 d 2 
ie — = - — if>(x,z) + U(x,z) i/)(x,z), (1) 

where z = ct is the longitudinal coordinate of the beam motion and U(x, z) = u(x, z)/m , yf3 2 c 2 
is the dimensionless potential. The transverse beam distribution function is given by 
a modulus squared of the beam wave function, p(x,z) = N \ip(x, z)\ 2 with a particle 
number of N. The beam wave function satisfies the following normalization condition, 
S-oo \ip(x,z)\ 2 dx = 1. 

The initial Gaussian beam wave function can be described by 



V*M)=(-y exp 
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where e = 2e rms with the unnormalized rms emittance of e rms . The parameters, <j\ and p\ are 
related to the twiss parameters of the initial particle distribution through o\ = y/e vms (3\ and 
Pi = —f3\/ot\ [5]. An elegant way to solve the differential equation (Eq. ([I])) is a functional 
integral method where the solution is given by the product of kernel (or propagator) and 
the initial beam wave function, 

/oo 
dx 1 K(x 2 , z; x u 0) ip(xx, 0). (3) 
-oo 

When we include the sextupole and octupole perturbation, the potential can be de- 
scribed as 



V(x) = U (x)-U p (x), (4) 

with 

U (x) = ^x 2 , U p (x) = ^ + ^x\ (5) 

where U (x) is the potential of a quadrupole magnet and U p (x) describes the aberration 
effects. The k 2 , k 3 , and fc 4 are related to the field strength of the quadrupole, sextupole, 
and octupole magnets. 



In the first-order perturbation, the kernels for a focusing and defocusing quadruple 
magnets become 



K f (x 2 , z;xi, 0;V p (x)) = Kf\x 2 , z;xi, 0) ( 1 + ^ h n (x 2 , z) 

\ n=0 

K d (x 2 ,z;x 1 ,0;V p (x)) = Kf\x 2 , z;x u 0) 1 + 9n{x 2 , z) 
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where Kf \x 2 , z;x±, 0) is the unperturbed kernel for a focusing case as given in Ref. [5]: 



Kf\x 2 ,z;xi,0) 
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The kernel for the defocusing case is obtained by replacing cot and esc functions in the 
focusing case with coth and csch functions, respectively. From the functional perturbation 
theory [7] , we obtained the coefficient functions for a focusing quadrupole as follows, 
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where r 3 = /c 3 <7o/(6/ci) and r 4 = fc40"g/(24/i;i) with ctq = e/{2y/k\). As noted in Ref. pE], r 3 
and r4 are the perturbation parameters. The functions Si(z) are given by 

si(z) 
s 2 (z) 

s 4 (;z) 

S5(«) 

s 7 (z) = 12y^^ - 8sin(2^z) +sin(4y^z). (10) 

The corresponding functions 2) for a defocusing case can be obtained by replacing the 
trigonometric functions in ht(x, z) with the corresponding hyperbolic functions. For a drift 
space, the unpertubed kernel is enough to get the final wave function: 



tan 2 (^£4,2/2), 

sec 2 (\fkiz/2) taxi(\Jkiz / 2) , 

(2 + sec 2 ( \fk[z/2)) txa.(<Jk x z/2), 

csc(Jkiz), 

k\[k~iz + 2\J~k[z cos{2\Jkiz) — 3sm.{2\Jkiz), 
\2\jk\z cos(\Jhi_z) — 9 sin(J~kiz) — sin(3-y/fciz) 
\2Jk x z - 8 sin(2\/%z) + sia(AJkxz). 
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where z is the length of the drift space. 

Because the beam wave function includes some perturbed terms after passing a quadrupole 
magnet with aberration, we generalize the initial beam wave function like 
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where a n = 0(k 3 , fc 4 ) describes the perturbation effects. 
After a drift space, the wave function becomes 



(IxiKq(x, z; xi, 0)ipi(xi, 0) 
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where the parameters 82, cr 2 , and P2 after the drift space are related to the corresponding 
initial parameters 9i,a±, and p±. The explicit formula can be found in Ref. |5j. The 
functions I n (x, 9) in the final wave function are given by 



Io(x, 



1, I 1 (x,9) = h™H 1 (x), I 2 ^e) = - A e m H 2 {x)+ 1 -, 
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where H n (x) are the Hermite polynomials. 

After a focusing lens of the effective length of z, the wave function becomes 
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For a defocusing case, the function h n is replaced with g n . The initial and final parameters, 
9i, <Ti, and pi, can be obtained by the same formula as given in Ref. [5]. 

We note that the relation between initial model parameters, o\ and p±, and result- 
ing parameters, 02 and P2, can be obtained by a well-known transformation of the twiss 
parameters [8j, 



(16) 



with jiPi — a\ = 1. The twiss parameters are related to the TWM paramters through 
= V^rmsPi and pi = -AM- 

We also note that the final wave function is not normalized to one even though the initial 
wave function is normalized to be one. The reason is that the modulus squared of the wave 
function includes the second-order terms which are proportional to 0{k%, /c|). Because the 
particle distribution function becomes negative in some region if we keep the linear terms 
only in modulus squared, we have to use a modulus squared for the particle distribution. 
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Hence the normalization increase should be an indicator to evaluate the validity of the 
perturbation calculation. We found that normalization increases less than few percent 
from one, the perturbation calculation described the PARMILA calculation successfully. 

In order to check the perturbation calculation, we compared the results with the PARMILA 
simulation with 50,000 macro particles which pass through a FODO lattice. The field gradi- 
ent and the effective length of a quadrupole magnet are 10 T/m and 0.2 m, respectively. The 
length of a drift space is 0.5 m. In the simulation, we assumed that the field strengths of the 
sextupole and octupole components are 3% of the quadrupole strength. We used matched 
input beam of the lattice with the twiss parameters, a = —1.60, f3 = 2.37 mm/mrad, and 
e r ms = 1-56 7r mm-mrad. They correspond to U\ — 1.92 mm and p 1 = 1.48 m. We used a 
Gaussian input beam and neglected the space charge effects. 

Figure 1 shows the input particle distribution (histogram) for the PARMILA calculation 
and the model result (solid line). The input beam wave function is normalize to one in the 
mks unit. The particle distributions after each beam-optical component are given in from 
Figure 2 to Figure 5. The solid lines in figures show the TWM result with sextupole and 
octupole perturbations. We also included the pure quadrupole result as dotted lines in 
the figures to evaluate the validity of the results with and without aberration effects. The 
PARMILA calculation shows clearly that particle distribution deviates from a Gaussian 
distribution. The TWM with aberration successfully described the asymmetric result as 
shown in the figures. In order to compare the results quantitatively, we calculated x values 
in each step. The result is summarized in Table I. It shows that the TWM with aberration 
improves the result than the model without perturbation. 

This work is related to the TWM calculation for a transverse beam dynamics of the 
charged particle with a sextupole and octupole perturbations. This model calculation is 
valid if space charge effects are negligible. Even though the perturbation result is lim- 
ited to the first-order, it explained the PARMILA simulation results successfully when the 
aberration field strength is less than few percent of the quadrupole strength. 
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Table 1: x values between PARMILA result and TWM calucation with and without aber- 
ration effects. 
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Figure 1: Particle distributions of the input beam. The histogram and solid line represent 
the PARMILA results and the model predictions, respectively. 
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Figure 2: Particle distributions after the first quadrupole magnet. The histogram, solid, and 
dotted lines represent the PARMILA results and the model predictions with and without 
aberration, respectively. 
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Figure 3: Particle distributions after the first drift space. The histogram, solid, and dotted 
lines represent the PARMILA results and the model predictions with and without aberra- 
tion, respectively. 
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Figure 4: Particle distributions after the second quadrupole magnet. The histogram, solid, 
and dotted lines represent the PARMILA results and the model predictions with and with- 
out aberration, respectively. 




Figure 5: Particle distributions after the second drift space. The histogram, solid, and 
dotted lines represent the PARMILA results and the model predictions with and without 
aberration, respectively. 
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